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Abstract We investigate the superfluid phase transition and effects of mass imbal- 
ance in the BCS (Bardeen-Cooper-Schrieffer)-BEC (Bose-Einstein condensation) 
crossover regime of an cold Fermi gas. We point out that the Gaussian fluctua- 



tion theory developed by Nozieres and Schmitt-Rink and the T-matrix theory, that 
are now widely used to study strong-coupling physics of cold Fermi gases, give 
unphysical results in the presence of mass imbalance. To overcome this problem, 
we extend the T-matrix theory to include higher-order pairing fluctuations. Using 
this, we examine how the mass imbalance affects the superfluid phase transition. 
Since the mass imbalance is an important key in various Fermi superfluids, such 
as 40 K- 6 Li Fermi gas mixture, exciton condensate, and color superconductivity in 
a dense quark matter, our results would be useful for the study of these recently 
developing superfluid systems. 

PACS numbers: 74.70.Tx,74.25.Ha,75.20.Hr 



1 Introduction 

The realization of a superfluid Fermi gas with mass imbalance is an exciting chal- 
lenge in cold atom physics. Once this pairing state is achieved, using a tunable 
interaction associated with a Feshbach resonance, we can study various properties 
of this novel superfluid phase from the weak-coupling BCS regime to the strong- 
coupling BEC limit in a unified manner 6 . Indeed, the so-called Sarma phase 4 
discussed in a spin-polarized Fermi superfluid has also been predicted in this sys- 
tem 5 . Since the formation of Cooper pairs between different species has been 
discussed in various systems, such as an exciton condensate in a semiconduc- 
tor^, and color superconductivity in a dense quark matter 3 , the realization of 
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mass imbalanced superfluid Fermi gas would contribute to the further develop- 
ments of these recently developing research fields. So far, the superfluid phase 
has not been realized in a mass-imbalanced Fermi gas. However, the formation 
of hetero-molecules was recently achieved in a 40 K- 6 Li Fermi-Fermi mixture-^, so 
that their superfluid phase transition might be realized near future. 

In this paper, we investigate the superfluid phase transition of an ultracold 
Fermi gas with mass imbalance in the BCS-BEC crossover region. In the ordi- 
nary mass-balanced case, the strong-coupling theory developed by Nozieres and 
Schmitt-Rink (NSR) 9 , as well as the T-matrix approximation (TMA) 10 ' 11 ' 12 , have 
been extensively used to successfully explain various aspects of the BCS-BEC 
crossover physics. However, we show that these strong-coupling theories do not 
work well in the presence of mass imbalance. A similar problem has also been 
pointed out in a superfluid Fermi gas with population imbalance 161 - . In this pa- 
per, we extend the T-matrix theory to include higher order pairing fluctuations, 
so as to overcome this serious problem. Using this extended T-matrix theory, we 
examine how the mass imbalance affects the superfluid phase transition in the 
BCS-BEC crossover region. 



2 Extended T-matrix for a Fermi gas with mass imbalance 

We consider a uniform two-component Fermi gas, described by the Hamiltonian, 

H = L %P° C pa c po - L C l+q/2t C -p+q/2l C -p'+l/n C p'+l/2- W 

P,o 1 p,p' 

Here, c po is an annihilation operator of a Fermi atom with momentum p and 
pseudospin a =t,i describing two atomic hyperfine states. E, pa = p 2 / {2m a ) — jx a 
is the kinetic energy of the a-spin component, measured from the Fermi chemical 
potential (where m a is the mass of a a-spin atom). The pairing interaction —U 
(< 0) is assumed to be tunable by adjusting the threshold energy of a Feshbach 
resonance (although we have ignored the detailed Feshbach mechanism in Eq. 
([T]l). As usual, we measure the interaction strength in term of the inverse scattering 
length a s , defined by 

4na, -U 
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where m is twice the reduced mass (2m~ l = m^ 1 +m^ 1 ). In this scale, the weak- 
coupling BCS regime and the strong-coupling BEC regime are, respectively, char- 
acterized by (&f0j) _1 — 1 and 1 > (^pc^) -1 (where k-p = (37T./V) 1 / 3 is the Fermi 
momentum of a gas of N Fermi atoms). The region between — 1 <^ (k^a s )~ l < 1 
is called the crossover region. 

To overcome the above mentioned breakdown of the NSR theory and T-matrix 
theory in the presence of mass imbalance, we extend the T-matrix theory to in- 
clude higher order fluctuations in the Cooper channel. In our extended T-matrix 
approximation (ETMA), the self-energy ,E CT (p,/oty,) in the single-particle Green's 
function, 

Gpa = ~. ? / ; Tj (3) 
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Fig. 1 (a) particle-particle scattering matrix F(q,i'v„) in the ladder approximation, (b) self- 
energy L a {p, ico„) in the extended 7-matrix theory, (c) self-energy E™ A (p, ia>„) in the ordinary 
T-matrix theory. In this figure, the double solid line and the single solid line describe G and G°, 
respectively. The wavy line denotes the pairing interaction —U. 



has the form 

E a (p,i(O n ) = T £ r(q,iv n )G q - p - a (iV„ - i(O n ). (4) 

Here, (O n and V„ are the fermion and boson Matsubara frequencies, respectively. 
The particle-particle scattering matrix, 

r ^ = - i-un iq , iVn y (5) 

involves effects of pairing fluctuations within the ladder approximation, which is 
diagrammatically described by FigQIa). Here, G° pa {i(O n ) = [i(O n — ^po-] -1 is the 
Green's function for a free Fermi gas, and 

n(q, iv n ) = T £ G° p+q/2 it (fv„ + i(0„)G°_ p+q/u (-i(O n ) (6) 
pja>„ 

is the lowest-order pair-correlation function. 

In the ordinary (non-self-consistent) T-matrix approximation (TM A), the Green's 
function G in Eq.Q is replaced by the non-interacting Green's function as 

E™ A (p,ico n ) = T £ r(q,iv n )G Q q _ p ,_ a (iv n - ion). (7) 

q.v„ 

Equation (0) is also used in the NSR theory, but the Green's function in Eq. (O is 
expanded up to 0(E™ A ) as 

G^ R (/ft)„) = G° po (i(O n ) + G%{i(O n )I™ A {p, i(o n )G%(i(O n ). (8) 

The superfluid phase transition temperature T c is conveniently determined from 
the Thouless criterion, which states that the phase transition occurs when the 
particle-particle scattering matrix F(q, iv n ) has a pole at q = V„ = 0. Since the 
three approximate theories, ETMA, TMA, and NSR, all use the same F(q, iv n ) in 
Eq. ©, the resulting r c -equation is also the same among them, which is given by 

= U y tanh(^ pt /(2r))+tanh(^ pi /(2r)) 

2 V Srt + Spi 
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Fig. 2 (Color online) Calculated superfluid phase transition temperature T c in the BCS-BEC 
crossover, normalized by Tp = (37r 2 A'/2v / 2) 2/ ' 3 m~ l . (a) ordinary T-matrix approximation 
(TMA). (b) extended T-matrix approximation (ETMA). The ratio of the mass imbalance is taken 
to be m-f /m± = 0.9. 



Since the chemical potential /J. a is known to remarkably deviate from the Fermi 
energy in the BCS-BEC crossover, we actually solve the r c -equation (|9), together 
with the equations for the number N a= + j_(= N/2) of a-spin atoms, 

j = T £ G pa (i(0 n ), (10) 

pd<o„ 

and determine T c , fU, jUj,, self-consistently. 



3 Superfluid phase transition ans effects of mass imbalance 

Figure[2a) shows the superfluid phase transition temperature T c for m-j-/mj = 0.9, 
calculated within TMA (where the self-energy in Eq. (0 is used). Clearly, the 
calculated T c exhibits unphysical behavior around the unitarity limit ((fcp<2.s) _1 ~ 
0). A similar unphysical result for T c is also obtained in the NSR theory (although 
we do not explicitly show the result here). 

To understand the origin of this unphysical result, it is helpful to compare 
the present system with a polarized Fermi gas, where the breakdown of the NSR 
theory and TMA have also been pointed ou t 15 i 16 i 14 . For this purpose, we rewrite 
the kinetic energies ff= -M as, 

p - ul( eL 

5pt m t \2m 

p = — (eL 

p ^ mi \2m 

where ji = (ju-j- + jUj.)/2, and 



■HJ-h, 



(11) 
(12) 



(13) 
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Fig. 3 (Color online) Calculated ETMA chemical potential [i„ in the BCS-BEC crossover at T c , 
normalized by the Fermi energy £p = (3n 2 N /2\/2) 2/ ' 3 m~ 1 of a Fermi gas with mass m. We set 
m-f/m^ = 0.9. 

When h = 0, the Fermi surfaces of the two components, determined by ^ p CT = 0, 
coincide with each other. On the other hand, the mismatch of the Fermi surfaces 
occurs when h 7^ 0. That is, h effectively works as an external magnetic field ap- 
plied to the system. In this sense, effects of mass imbalance is similar to those 
of population imbalance (where the mismatch of the Fermi surfaces between the 
two components is induced by setting AT+ 7^ A^.) In the latter spin-polarized sys- 
tem, the reason for the breakdown of the NSR theory and TMA is known to be 
the incomplete treatments of magnetic fluctuations and pseudogap effects associ- 
ated with pairing fluctuation s 15 i 16 i 14 . It has also been pointed out 14 in the case 
of spin-polarized Fermi gas that one can eliminate this serious problem when 
the self-energy in FigQJb) is used. Then, because of the similarity between the 
two systems, ETMA is expected to be also valid for the problem in the mass- 
imbalanced case. Indeed, as shown in Fig|2b), ETMA gives the expected smooth 
crossover behavior of T c in the BCS-BEC crossover. That is, starting from the 
weak-coupling regime, T c gradually increases with the increase of the interaction 
strength, and it approaches a constant value in the BEC regime. In the BEC limit, 
the system is well described by a gas of tightly bound hetero-molecules with the 
particle number N/2 and molecular mass M = m^ + m±. Noting this, one finds 

2tn 

T c = ■ x0.2187f (eeTbec), (14) 

'"• • '". 

in the BEC limit (where Tf = (3n 2 N /2\/l) 2 l 3 m~ x is the Fermi temperature of a 
Fermi gas with the mass m). 

Figure [3] shows the ETMA chemical potential fj. a in the BCS-BEC crossover 
at T c . As expected from the ETMA result on T c in Fig|2tb), one sees smooth be- 
havior of the chemical potential as a function of the interaction strength. As in the 
mass-balanced case, both /J| and /i| gradually decrease, as one passes through the 
BCS-BEC crossover region. They become negative in the BEC regime, which is a 
typical phenomenon of the BCS-BEC crossover physics. 

We now examine effects of the mass imbalance on the superfluid phase transi- 
tion temperature T c by using the extended T-matrix theory. As shown in Fig|4] 
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Fig. 4 (Color online) Superfluid phase transition temperature T c as a function of the mass imbal- 
ance ratio m^/my. (a) extended J-matrix approximation (ETMA). (b) mean-field theory (MF). 
In this figure, Tbec is the BEC phase transition temperature of an ideal hetero-molecular Bose 
gas, given by Eq. dl4l . 



the mass imbalance tends to decrease T c , irrespective of the magnitude of the 
interaction strength. In the BEC regime ((&f«s) 1 = 2), a finite T c is obtained 
down to the limit m+ jm^ — > 0. In this regime, hetero-molecules have already been 
formed above T c , the superfluid phase transition is dominated by the BEC of these 
molecules. Indeed, as shown in FigHa), T c in this regime is well described by the 
BEC phase transition temperature of an ideal Bose gas, given by Tbec in Eq. (114b . 

In contrast to the BEC regime, T c vanishes at a certain value of m^/m^ (< 1) 
in the BCS regime, as shown in Fig.|4ja). At a glance, this vanishing T c seems to 
be because of the depairing effect by the effective magnetic field h in Eq. < l 1 3 b - 
However, when we substitute the Fermi energy E^o = (3n 2 N /2v / 2) 2 ' 3 wi^ 1 of the 
<7-spin component into pL a in Eq. ( 1131) . we find h = 0. This is simply because 
the magnitude of the Fermi momentum in a free Fermi gas is independent of the 
particle mass at T = 0. Actually, h becomes finite at finite temperatures, because 
the temperature dependence of the chemical potential usually depends on m a . 
However, when we ignore all many-body fluctuation effects and simply treat the 
temperature effect within the mean-field theory, a finite T c is obtained everywhere 
irrespective of the magnitude of m+/m^, as shown in Figj4jb). This is simply be- 
cause the mismatch of the Fermi surfaces between the two components always 
disappears at T = within the mean-field level. Thus, the vanishing T c seen in 
Figffia) ls considered to come from an interaction effect beyond the mean-field 
level. Indeed, when the Fermi surface is smeared by a many-body effect, it is ex- 
pected to give an effect similar to the temperature, leading to a finite h. 

We briefly note that, while the present ETMA involves higher order fluctuation 
effects than the ordinary T-matrix theory, it is still not fully self-consistent in the 
sense that the non-interacting Green's function G° is used in the pair correlation 
function in Eq. ©. In addition, we have assumed the second order phase transi- 
tion, as well as the simplest uniform s-wave pairing state, in determining T c . Thus, 
since the vanishing T c , as well as the first order phase transition like behavior of T c 
seen in Fig|4jb), might depend on details of the theory, further analyses on these 
would be necessary, which remains as our future problem. 
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4 Summary 

To summarize, we have investigated the superfluid phase transition of a cold Fermi 
gas with mass imbalance. We have extended the T-matrix theory to include higher 
order fluctuation effects, so as to overcome the serious problems existing in the 
NSR theory and the T-matrix theory. Using this extended T-matrix theory, we 
have calculated T c in the BCS-BEC crossover region to examine effects of mass 
imbalance. We show that, within ETMA, while a finite T c is always obtained in the 
BEC regime, T c vanishes at a certain value of m^/m^ in the weak-coupling BCS 
regime. Since the Fermi superfluids/condensates have been discussed, not only in 
the field of cold Fermi gas physics, but also in semiconductor physics (exciton), 
as well as in high-energy physics (color superconductivity), our results would be 
useful for the study of these hetero-type Fermi superfiuids. 
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